New-generation electron-propagator methods for vertical electron detachment energies of molecular anions: benchmarks and applications to model green-fluorescent-protein chromophores.

Ab initio electron-propagator calculations continue to be useful companions to experimental investigations of electronic structure in molecular anions. A new generation of electron-propagator methods recently has surpassed its antecedents' predictive accuracy and computational efficiency. Interpretive clarity has been conserved, for no adjustable parameters have been introduced in the preparation of molecular orbitals or in the formulation of approximate self-energies. These methods have employed the diagonal self-energy approximation wherein each Dyson orbital equals a canonical Hartree-Fock orbital times the square root of a probability factor. Numerical tests indicate that explicitly renormalized, diagonal self-energies are needed when Dyson orbitals have large valence nitrogen, oxygen or fluorine components. They also demonstrate that even greater accuracy can be realized with generalizations that do not employ the diagonal self-energy approximation in the canonical Hartree-Fock basis. Whereas the diagonal methods have fifth-power arithmetic scaling factors, the non-diagonal generalizations introduce only non-iterative sixth-power contractions. Composite models conserve the accuracy of the most demanding combinations of self-energy approximations and flexible basis sets with drastically reduced computational effort. Composite-model results on anions that resemble the chromophore of the green fluorescent protein illustrate the interpretive capabilities of explicitly renormalized self-energies. Accurate predictions on the lowest vertical electron detachment energy of each anion confirm experimental data and the utility of the diagonal self-energy approximation.

New-Generation Electron-Propagator Methods for Molecular Electron-Binding Energies.

A new generation of electron-propagator methods for the calculation of electron binding energies has surpassed its antecedents with respect to accuracy, efficiency, and interpretability. No adjustable parameters are introduced in these fully ab initio procedures. Numerical tests versus several databases of valence, vertical electron binding energies of closed-shell molecules and atoms have been performed. Easily interpreted self-energy approximations with cubic arithmetic scaling produce mean absolute errors (MAEs) of 0.2 and 0.3 eV for electron detachments and attachments, respectively. The most accurate explicitly renormalized methods with fifth-power arithmetic scaling yield MAEs below 0.1 eV for detachments and attachments. Approximate renormalization leads to more efficient fifth-power alternatives for electron detachments that achieve similar accuracy with fewer bottleneck operations. Composite protocols generate excellent predictions versus highly accurate basis-extrapolated standards and experiments. The validity of the diagonal self-energy approximation and the accuracy of the approximate renormalizations are confirmed. The success of these perturbative methods based on canonical Hartree-Fock orbitals rests on a Hermitized, intermediately normalized superoperator metric. The results of all of the new-generation calculations may be analyzed in terms of final-state orbital relaxation and differential correlation effects.

New-Generation Electron-Propagator Methods for Calculations of Electron Affinities and Ionization Energies: Tests on Organic Photovoltaic Molecules.

A new generation of ab initio electron-propagator self-energies recently superseded its antecedents' accuracy and computational efficiency in calculating vertical ionization energies (VIEs) of closed-shell molecules. (See J. Chem. Phys. 2021, 155, 204107, J. Chem. Theory Comput. 2022, 18, 4927, J. Chem. Phys. 2023, 159, 124109.) No adjustable parameters were introduced in the generation of reference orbitals or in the construction of self-energies. The same approach has been extended in this work to vertical electron affinities (VEAs). Calculations were performed on 24 conjugated, organic photovoltaic molecules with diverse functional groups. These molecules are considerably larger than those studied in previous tests on VIEs. Several new-generation self-energies produce mean absolute errors (MAEs) below 0.1 eV versus ΔCCSD(T) (i.e., total energy differences from the coupled-cluster singles, doubles, and perturbative triples method) VIEs and VEAs obtained with identical basis sets. A composite model employs cubically and quintically scaling algorithms and power-law basis-set extrapolations based on augmented double-triple or triple-quadruple ζ data. Its MAEs are near 0.05 eV versus benchmark values, with 0.03 eV error bars for the lowest VIE and the highest VEA of each molecule. A more efficient and equally accurate composite model for calculating VIEs avoids full transformations of electron repulsion integrals to the molecular orbital basis. High probability factors support the diagonal self-energy approximation, wherein Dyson orbitals are proportional to canonical, Hartree-Fock orbitals.

A new generation of non-diagonal, renormalized self-energies for calculation of electron removal energies.

A new generation of diagonal self-energies for the calculation of electron removal energies of molecules and molecular ions that has superseded its predecessors with respect to accuracy, efficiency, and interpretability is extended to include non-diagonal self-energies that permit Dyson orbitals to be expressed as linear combinations of canonical Hartree-Fock orbitals. In addition, an improved algorithm for renormalized methods eliminates the convergence difficulties encountered in the first studies of the new, diagonal self-energies. A dataset of outer-valence, vertical ionization energies with almost full-configuration-interaction quality serves as a standard of comparison in numerical tests. The new non-diagonal, renormalized methods are slightly more accurate than their diagonal counterparts, with mean absolute errors between 0.10 and 0.06 eV for outer-valence final states. This advantage is procured at the cost of an increase in the scaling of arithmetic bottlenecks that accompany the inclusion of non-diagonal self-energy terms. The new, non-diagonal, renormalized self-energies are also more accurate and efficient than their non-diagonal predecessors.

Electron Propagator Theory of Vertical Electron Detachment Energies of Anions: Benchmarks and Applications to Nucleotides.

A new generation of ab initio electron-propagator self-energy approximations that are free of adjustable parameters is tested on a benchmark set of 55 vertical electron detachment energies of closed-shell anions. Comparisons with older self-energy approximations indicate that several new methods that make the diagonal self-energy approximation in the canonical Hartree-Fock orbital basis provide superior accuracy and computational efficiency. These methods and their acronyms, mean absolute errors (in eV), and arithmetic bottlenecks expressed in terms of occupied (O) and virtual (V) orbitals are the opposite-spin, non-Dyson, diagonal second-order method (os-nD-D2, 0.2, OV2), the approximately renormalized quasiparticle third-order method (Q3+, 0.15, O2V3) and the approximately renormalized, non-Dyson, linear, third-order method (nD-L3+, 0.1, OV4). The Brueckner doubles with triple field operators (BD-T1) nondiagonal electron-propagator method provides such close agreement with coupled-cluster single, double, and perturbative triple replacement total energy differences that it may be used as an alternative means of obtaining standard data. The new methods with diagonal self-energy matrices are the foundation of a composite procedure for estimating basis-set effects. This model produces accurate predictions and clear interpretations based on Dyson orbitals for the photoelectron spectra of the nucleotides found in DNA.

Electron Propagator Self-Energies versus Improved GW100 Vertical Ionization Energies.

Ab initio electron propagator (EP) methods that are free of adjustable parameters in their self-energy formulae and in the generation of their orbital bases have been applied to the calculation of the lowest vertical ionization energies (VIEs) of the GW100 set. An improved set of standard results accompanied by irreducible representation assignments has been produced indirectly with coupled-cluster singles and doubles plus perturbative triples, i.e., CCSD(T), total energy differences at initial-state geometries reoptimized (in 28 cases) with the largest applicable point groups. The best compromises of accuracy and efficiency belong to a new generation of EP self-energies, several members of which may be derived from an intermediately normalized, Hermitized super-operator metric. The following diagonal self-energy methods are optimal: opposite-spin non-Dyson second order (os-nD-D2), approximately renormalized partial third order (P3+), approximately renormalized quasiparticle third order (Q3+), and non-Dyson approximately renormalized linear third order version B (nD-L3+B). Their mean absolute errors (MAEs) in electron volts and arithmetic scaling factors expressed in terms of occupied (O) and virtual (V) orbital dimensions are, respectively, (0.18, OV2), (0.14, O2V3), (0.15, O2V3), and (0.11, OV4). The 0.06 eV MAE for the non-diagonal, sixth-power (O2V4) Brueckner doubles, triple-field operator (BD-T1) EP method is exceeded by the 0.1 eV MAE with respect to experiments in seventh-power, ΔCCSD(T) calculations and indicates that BD-T1 may serve as a direct, spin-symmetry-conserving alternative in the generation of standard results for VIEs of larger, closed-shell molecules.

A new generation of diagonal self-energies for the calculation of electron removal energies.

A new generation of diagonal self-energy approximations in ab initio electron propagator theory for the calculation of electron removal energies of molecules and molecular ions has been derived from an intermediately normalized, Hermitized super-operator metric. These methods and widely used antecedents such as the outer valence Green's function and the approximately renormalized partial third order method are tested with respect to a dataset of vertical ionization energies generated with a valence, triple-ζ, correlation-consistent basis set and a converged series of many-body calculations whose accuracy approaches that of full configuration interaction. Several modifications of the diagonal second-order self-energy, a version of G0W0 theory based on Tamm-Dancoff excitations and several non-diagonal self-energies are also included in the tests. All new methods employ canonical Hartree-Fock orbitals. No adjustable or empirical parameters appear. A hierarchy of methods with optimal accuracy for a given level of computational efficiency is established. Several widely used diagonal self-energy methods are rendered obsolete by the new hierarchy whose members, in order of increasing accuracy, are (1) the opposite-spin non-Dyson diagonal second-order or os-nD-D2, (2) the approximately renormalized third-order quasiparticle or Q3+, (3) the renormalized third-order quasiparticle or RQ3, (4) the approximately renormalized linear third-order or L3+, and (5) the renormalized linear third-order or RL3 self-energies.

Eigenvalues of uncorrelated, density-difference matrices and the interpretation of Δ-self-consistent-field calculations.

Two theorems on the eigenvalues of differences of idempotent matrices determine the natural occupation numbers and orbitals of electronic detachment, attachment, or excitation that pertain to transitions between wavefunctions that each consist of a single Slater determinant. They are also applicable to spin density matrices associated with Slater determinants. When the ranks of the matrices differ, unit eigenvalues occur. In addition, there are ±w pairs of eigenvalues where |w| ≤ 1, whose values are related to overlaps, t, between the corresponding orbitals of Amos and Hall, and Löwdin by the formula w=±1-t2 12. Generalized overlap amplitudes, including Dyson orbitals and their probability factors, may be inferred from these eigenvalues, which provide numerical criteria for: classifying transitions according to the number of holes and particles in final states with respect to initial states, identifying the most important effects of orbital relaxation produced by self-consistent fields, and the analysis of Fukui functions. Two similar theorems that apply to sums of idempotent matrices regenerate formulae for the natural orbitals and occupation numbers of an unrestricted Slater determinant that were published first by Amos and Hall.

Dyson-orbital concepts for description of electrons in molecules.

Dyson orbitals, their electron-binding energies, and probability factors provide descriptions of electrons in molecules that are experimentally verifiable and that generalize qualitatively useful concepts of uncorrelated, molecular-orbital theory to the exact limit of Schrödinger's time-independent equation. Dyson orbitals are defined as overlaps between initial, N-electron states and final states with N ± 1 electrons and therefore are useful in the prediction and interpretation of many kinds of spectroscopic and scattering experiments. They also are characteristic of N-electron initial states and may be used to construct electron densities, one-electron properties, and total energies with correlated Aufbau procedures that include probability factors between zero and unity. Relationships with natural orbitals, Kohn-Sham orbitals, and Hartree-Fock orbitals facilitate insights into the descriptive capabilities of Dyson orbitals. Electron-propagator approximations that employ the Dyson quasiparticle equation or super-operator secular equations enable direct determination of Dyson orbitals and obviate the need for many-electron wavefunctions of initial or final states. Numerical comparisons of the amplitudes and probability factors of Dyson orbitals calculated with several self-energy approximations reveal the effects of electron correlation on these uniquely defined, one-electron wavefunctions.

Electron Propagator Methods for Vertical Electron Detachment Energies of Anions: Benchmarks and Case Studies.

Ab initio electron propagator methods are efficient and accurate means of calculating vertical electron detachment energies of closed-shell, molecular anions with nuclei from the first three periods. Basis set extrapolations enable definitive comparisons between electron propagator results and benchmarks defined by total energy differences obtained with coupled-cluster, single, double, plus perturbative triple substitution theory. The best compromises of accuracy and efficiency are provided by the renormalized, partial third-order, diagonal (P3+) self-energy and by the nondiagonal, renormalized, second-order (NR2) approximation. The outer-valence Green function, the two-particle-one-hole Tamm-Dancoff approximation, the third-order algebraic diagrammatic construction, and the renormalized third-order methods also are examined. A detailed analysis of errors for small anions is performed. Case studies include F-(H2O) and Cl-(H2O) complexes, C5H5-, two P2N3- pentagonal rings, and a superhalide, Al(BO2)4-, whose electron detachment energy is more than double those of the halide anions. These applications illustrate the versatility of electron propagator methods, their utility for interpreting negative-ion photoelectron spectra, and their promise in the discovery of unusual properties and patterns of chemical bonding. Composite methods, which combine basis set effects calculated at the relatively efficient diagonal, second-order level and higher correlation effects calculated with small basis sets, provide excellent estimates of basis set-extrapolated P3+ or NR2 results and facilitate applications to large molecules. In the P3+ and NR2 methods, a judicious choice of low-order couplings between hole operators that correspond to the assumptions of Koopmans's theorem and operators that describe final-state relaxation and polarization and initial-state correlation leads to predictive accuracy, computational efficiency, and interpretive lucidity.

Communication: Explicitly correlated formalism for second-order single-particle Green's function.

We present an explicitly correlated formalism for the second-order single-particle Green's function method (GF2-F12) that does not assume the popular diagonal approximation and describes the energy dependence of the explicitly correlated terms. For small and medium organic molecules, the basis set errors of ionization potentials of GF2-F12 are radically improved relative to GF2: the performance of GF2-F12/aug-cc-pVDZ is better than that of GF2/aug-cc-pVQZ, at a significantly lower cost.

One-particle many-body Green's function theory: Algebraic recursive definitions, linked-diagram theorem, irreducible-diagram theorem, and general-order algorithms.

A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green's function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green's function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green's function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green's function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green's function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams.

Assessment of Electron Propagator Methods for the Simulation of Vibrationally Resolved Valence and Core Photoionization Spectra.

The analysis of photoelectron spectra is usually facilitated by quantum mechanical simulations. Because of the recent improvement of experimental techniques, the resolution of experimental spectra is rapidly increasing, and the inclusion of vibrational effects is usually mandatory to obtain a reliable reproduction of the spectra. With the aim of defining a robust computational protocol, a general time-independent formulation to compute different kinds of vibrationally resolved electronic spectra has been generalized to also support photoelectron spectroscopy. The electronic structure data underlying the simulation are computed using different electron propagator approaches. In addition to the more standard approaches, a new and robust implementation of the second-order self-energy approximation of the electron propagator based on a transition operator reference (TOEP2) is presented. To validate our implementation, a series of molecules has been used as test cases. The result of the simulations shows that, for ultraviolet photoionization spectra, the more accurate nondiagonal approaches are needed to obtain a reliable reproduction of vertical ionization energies but that diagonal approaches are sufficient for energy gradients and pole strengths. For X-ray photoelectron spectroscopy, the TOEP2 approach, besides being more efficient, is also the most accurate in the reproduction of both vertical ionization energies and vibrationally resolved bandshapes.

Comment on "Are polynuclear superhalogens without halogen atoms probable? A high-level ab initio case study on triple-bridged binuclear anions with cyanide ligands" [J. Chem. Phys. 140, 094301 (2014)].

For the vertical electron detachment energies of triply-bridged Mg2(CN)5- superhalides, the Outer Valence Green Function (OVGF) yields results similar to those of the coupled-cluster singles and doubles plus approximate triples, or CCSD(T), method. Invalid comparisons between states with different symmetry or localization properties underlie the assertion of Yin et al that OVGF produces large discrepancies with respect to CCSD(T) for several isomers of Mg2(CN)5-.

Composite electron propagator methods for calculating ionization energies.

Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green's function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C60 illustrate the capabilities of CEP methods for calculations on large molecules.

Comment on "Is the regulation of the electronic properties of organic molecules by polynuclear superhalogens more effective than that by mononuclear superhalogens? A high-level ab initio case study" by M.-M. Li, J.-F. Li, H.-C. Bai, Y.-Y. Sun, J.-L. Li and B. Yin, Phys. Chem. Chem. Phys., 2015, 17, 20338.

The Outer Valence Green Function (OVGF) and coupled-cluster singles and doubles plus approximate triples, or CCSD(T), methods yield similar results for the vertical detachment energies of superhalides studied recently by Li et al. The success of the OVGF method contradicts claims by Li et al. in their recent article.